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INTRODUCTION 


The  possibility  of  interactions  between  very-low-frequency  electro- 
magnetic waves  and  energetic  electrons  in  the  earth’s  magnetosphere  has 
lor, 2 been  recognized.  Dowden  (1962)  proposedthat  whistlers  could  be 
generated  through  a Doppler  - shifted  resonance  with  single  electrons  or 
electron  bunches  near  the  equator.  Brice  (1963)  outlined  a method  for  such 
an  interaction  with  a general  population  of  energetic  electrons,  since  the 
recuired  electron  packet  was  an  improbable  initial  condition.  Brice  also 
pointed  out  that  this  interaction,  a cyclotron  or  transverse  resonance, 
would  result  in  the  lowering  of  the  pitch-angle  of  the  electrons.  Hell  well 
(1967)  later  generalized  this  result  to  permit  explanation  of  rising  tones 
through  off-equatorial  interaction  regions. 

A number  of  observations  of  naturally  occurring  VLF  waves,  generally 
lightning-generated  w'histlers  or  plasmapause  hiss,  have  been  reported  as 
being  correlated  with  ele ctron  pre cipitation.  Rosenberg  et  al.  (1971)  and 
Foster  and  Rosenberg  (1976),  for  instance,  report  correlations  between 
VLF  waves  and  x-rays  observed  by  balloon-borne  detectors.  The  actual 
particles  have  also  been  observed  on  satellites.  Holzer  et  al.  (1974) 
presented  correlations  between  chorus  observed  near  the  outer  edge  of  the 
stable  trapping  region  and  precipitation  of  >45  keV  electrons  on  OGO  6. 

Attempts  have  been  made  to  study  propagation  and  amplification  of 
VLF  waves  through  generation  and  radiation  of  test  signals  from  various 
places.  A significant  effort  has  centered  around  the  Siple  Station,  Antartica, 
transmissions  and  receptions  at  Roberval,  Canada  (Helliwell  and  Katsvif r akis, 
1974).  These  efforts  have  led  to  significant  results,  but  up  to  now  no  one 
has  reported  direct  evidence  (i.  e.  , particles  or  x-rays)  of  actual  pitch-angle 
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scattering  of  the  electrons  in  an  interaction  with  the  waves.  A second  j 

j 

major  effort  of  this  type,  which  involved  VLF  transmissions  from  Alaska  i 

(see  Koons  et  al.  , 1976),  also  appeared  to  produce  no  observable  electron 
precipitation.  In  this  latter  project,  some  evidence  of  proton  precipitation 
correlated  with  transmissions  was  observed  (Koons,  1975),  although  it  is 
unclear  hav  protons  could  be  affected  by  the  waves  themselves.  Bell  (1976) 
has  proposed  generation  of  ULF  waves  from  VLF  waves  through  an  intermediate 
interaction  with  energetic  electrons.  The  periodicity  in  amplification  of  VLF 
signals  reported  by  Helliwell  and  Katsufrakis  (1974)  might  be  an  indication 
that  the  Bell  suggestion  is  valid.  No  ground  station  magnetic  field  data  were 
a'  e to  determine  whether  the  modulated  amplification  was  accompanied 

aves.  If  such  waves  are  generated,  they  might  be  responsible  for 
jn  precipitation,  if  it  occurs. 

Other  strong  coherent  sources  of  VLF  waves  are  available  and  are 
routinely  operating,  although  for  the  purposes  of  communication  and  naviga- 
tion. No  direct  observation  of  particle  precipitation  due  to  these  sources 
has  been  reported  previously,  although  there  is  some  evidence  of  enhanced 
E-region  ionization  resulting  from  high  power  HF  transmissions  (Wright, 

1975).  A possibility  advanced  by  Wright  was  the  establishment  of  a duct  with  i 

the  high  power  HF  transmission,  natural  generation  of  VLF  waves  within 
the  duct,  and  resultant  precipitation  of  energetic  electrons  to  produce  the 

enhanced  E-layer  ionization.  In  the  present  paper,  we  report  the  observa-  ^ 

' 

tion  of  recurrent  energetic  electron  precipitation  with  a number  of  unique  ] 

parameters  which  we  feel  positively  identifies  it  as  being  caused  by  the 
Russian  VLF  communications  station,  UMS. 
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INSTRUMENTATION 


The  data  presented  here  were  obtained  by  sets  of  180°  focussing 
magnetic  spectrometers  on  two  different  satellites  during  the  1968-1970 
time  period.  The  bulk  of  the  data  was  obtained  from  the  OVl-19  (1969-25C) 
magnetic  electron  spectrometers  which  have  been  discussed  previously 
(Vampola,  1971).  But  the  initial,  and  in  some  ways  the  best,  observations 
were  made  by  a pair  of  magnetic  spectrometers  flown  on  the  OVl-14 
(1968-26B)  during  April,  1968.  The  OVl-14  was  a U.  S.  Air  Force  Office 
of  Aerospace  Research  satellite  launched  into  a 100°  inclination  elliptical 
orbit  on  April  6,  1968.  Apogee  and  perigee  were  9970  km  and  550  km, 
respectively.  The  vehicle  was  spin- stabilized  at  about  10  RPM  u-ith  the  spin 
also  being  used  to  provide  a scanning  function  for  the  scientific  payload.  An 
on-board  magnetometer  was  used  to  obtain  magnetic  aspect  for  the  particle 
instrumentation.  Several  days  after  launch,  a failure  in  the  spacecraft  power 
system  battery  charging  circuit  caused  overcharging  and  rupture  of  the 
battery,  which  resulted  in  premature  termination  of  the  mission.  An  on- 
board tape  recorder  stored  up  to  four  hours  of  data;  seven  tape  dumps  were 
achieved  over  a four-day  period  prior  to  termination  of  the  mission. 

The  magnetic  spectrometers  on  the  OVl-14  were  similar  to  others 
previously  (and  since)  flown  on  other  Air  Force  satellites.  One  was  identi- 
cal to  the  low  energy  unit  on  OVl-19  except  for  the  magnetic  field  (380  guass 
on  the  OVl-14)  and  the  acceptance  angle  of  the  aperture.  The  other  unit 
was  the  OV3-3  ( 1966-70A ) backup  instrument  (Vampola,  1969)  with  a modified 
collimator.  The  field  in  this  instrument  was  1790  gauss.  Each  instrument 
contained  a background  detector  to  measure  bremsstrahlung  and  penetrating 
particle  responses.  Table  1 gives  the  center  energies  and  geometric-energy 
factors  for  the  channels  in  both  units. 
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TABLE  1 


Low  Ener 

gy  Spectrometer 

High  Energy  Spectrometer 

Channel 

Center  Energy 
(MeV) 

Geometric-Ene  rgy 
Factor 

( cm^  - s ter  -ke  V) 

Center  Ene  rgy 
(MeV) 

Geometric  - Ener  gy 
F acto  r 

(cm^-  ster-keV ) 

1 

. 035 

1.42 

. 531 

15.  8 

-> 

. 062 

1, 60 

. 809 

14.  9 

3 

. 0^5 

1. 70 

1.  10 

13.  7 

4 

. 132 

1.75 

1. 40 

12.  6 

5 

. 174 

1.  72 

1.  70 

11.  6 

6 

. 219 

1.63 

2.  00 

10.  3 

7 

. 267 

1.  53 

2.  30 

9.  18 

8 

. 318 

1.43 

2.  60 

8.  29 

1 
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OBSERVATIONS 

The  time  period  of  the  OVl-14  data  was  the  recovery  period  following 
a moderate  magnetic  storm  on  April  5,  1968.  Highly  monoenergetic  precipi- 
tation spikes  were  seen  in  the  electron  distributions  over  a narrow  range  of 
L values  for  all  low  altitude  passes  over  central  and  western  Asia.  Examples 
of  these  spikes  are  shown  in  Figure  I.  The  highly  monoenergetic  nature  of 
these  distributions  suggested  a resonant  interaction  with  a monochromatic 
wave.  The  total  precipitation  event,  in  each  instance,  extended  over  about 
. 15  L,  but  the  observable  flux  enhancements  at  a given  energy  were  much 
narrower.  Fig.  2 illustrates  the  relationship  of  precipitated  flux  versus  L 
at  fixed  energies  and  also  the  dispersion  in  L as  a function  of  energy  for  the 
flux  peak.  The  relatively  high  spin  rate  coupled  . eh  a long  time  constant 
for  count  rate  meters  at  low  count  rates  serious'  / degraded  the  accuracy 
with  which  pitch- angle  distributions  could  be  determined.  The  limited 
amount  of  data  available  from  the  OVl-14  satellite  precluded  a definitive 
study  of  the  mechanism  involved  in  this  phenomenon  although  our  original 
conjectures  appear  to  be  correct. 

Recently,  a large  amount  of  data  from  the  OVl-19  became  available 
for  analysis.  Because  of  a low  efficiency  factor  in  the  lowest  energy  channel 
and  an  early  in-flight  failure  of  the  detector  in  the  next  lowest  energy  channel, 
the  OVl-19  data  set  is  not  ideal  for  a complete  study  of  the  energy  dependence 
of  these  precipitations.  However,  since  it  is  an  extensive  data  set,  a search 
was  made  for  events  of  the  type  seen  in  the  OVl-14  data  set.  Figure  3 is 
a typical  example  of  the  same  phenomenon  as  observed  by  the  OVl'  19  instru- 
mentation. Note  that  the  same  characteristics  are  present:  narrow  distribu- 
tion in  L at  any  given  energy;  decreasing  L location  with  increasing  energy: 

L value  between  about  1.6  and  1.85. 
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DAY 

97 

99 

99 


U.T. 

B 

L 

ALT 

E.L. 

LT. 

62,111 

0.395 

1.69 

714 

118.4“ 

1;03 

62,983 

0.299 

1.68 

1210 

127.5“ 

2:00 

63,945 

0.387 

1.68 

723 

109.6“ 

1:04 

Energy  spectra  obtained  by  magnetic  electron  spectrometers 
on  the  OVl-14  satellite  in  1968  in  the  drift  loss-cone  near 
L = 1.68.  The  narrow  energy  distributions  indicate  a resonant 
interaction  is  responsible  for  their  appearance  in  the  loss-cone 


:ure  2 Flux  vs.  time  plots  of  dcta  similar  to  those  which  produced  the 
narrow  energy  distributions  of  Figure  1.  Note  the  monotonic 
variation  in  L of  the  p<‘ ai<  of  the  precipitation  distribution  as  a 
function  of  energy.  The  scatter  in  the  data  points  is  the  result 
of  pitch-angle  sampling. 
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Figure  3 Data  similar  to  those  of  Figure  2,  but  obtained  by  the  OVl-19 
satellite  in  1969. 
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For  the  purposes  of  this  investigation,  use  was  made  of  the  fact 
that  electrons  with  small  equatorial  pitch- angles  will  be  lost  into  the 


earth's  atmosphere,  as  they  drift  eastward  due  to  geomagnetic  gradient 
effects,  at  a location  determined  by  their  pitch-angle  and  the  local  magneto- 
spheric  topography.  .A  thorough  discussion  of  this  effect  is  available 
elsewhere  (Roederer,  1970),  but  a brief  discussion  will  be  presented  here 
since  it  is  central  to  the  identification  of  the  source  of  these  precipitation 
events.  Figure  4 depicts  the  relationship  between  the  equatorial  pitch- 
angle  of  a particle  on  the  L = 1.75  shell  and  its  destiny  as  a function  of 
east  longitude.  For  the  purposes  of  Figure  4,  it  is  assumed  that  a particle 
that  has  a local  mirror  point  at  or  below  100  km  will  be  immediately  lost 
into  the  atmosphere;  i.  e.  , it  is  in  the  bounce  loss-cone.  A particle  which 
has  an  eauatorial  pitch-angle  between  that  of  the  local  bounce  loss-cone 
and  the  maxim.um  bounce  loss -cone  (which  in  this  case  occurs  at  about 
330°  E.  L.  ) is  defined  as  being  in  the  drift  loss-cone,  since  at  some  point 
in  its  drift,  it  will  encounter  a 100  km  mirror  altitude  and  be  lost  into  the 
atmosphere.  Particles  with  equatorial  pitch-angles  greater  than  this 
maximum  bounce  loss-cone  angle  are  considered  stably  trapped.  By 
limiting  observations  to  that  particle  population  which  is  currently  in  the 
drift  loss  cone,  one  can  identify  particles  that  have  been  pitch-angle 
scattered  during  the  immediately  preceeding  drift  period  (or  a portion  of 
it).  To  put  it  picturesquely,  for  the  example  in  Figure  4,  the  geomagnetic 
configuration  at  330°  E.  L.  has  wiped  the  slate  clean  for  equatorial  pitch- 
angles  up  to  27.6°.  Any  perturbation  at  any  longitude  which  lowers  the 
pitch-angle  of  a particle  to  below  27.6°  will  cause  the  particle  to  either 
be  locally  precipitated  or  to  reside  in  the  drift  loss-cone  until  it  encounters 
the  atmosphere  in  its  eastward  drift. 
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EQUATORIAL  PITCH  ANGLE 
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The  portion  of  the  OVl-19  data  set  that  was  available  for  this  investi- 
gation was  examined  for  the  presence  of  events  similar  to  those  of  Figs.  2 


and  3,  that  is,  of  particles  in  the  drift  loss  cone  for  the  L interval  1.  5 to 
2.0.  Of  the  eighty-eight  periods  of  observation  in  this  drift  loss-cone 
region,  thirty  precipitation  spikes  were  observed.  The  data  set  is  shown 
graphically  in  Fig.  5 with  time  and  east  longitude  as  the  parameters. 

Three  symbols  are  used  to  depict:  a)  observation  of  a precipitation  spike; 

b)  no  spike  present  although  the  spacecraft  was  in  a locatio n where  others 
had  been  observed;  c)  no  spike  present  but  the  spacecraft  was  below  the 
drift  shell  defined  by  the  atmosphere  at  50°  east  longitude.  (The  actual 
minimum  equatorial  pitch-angle  of  this  drift  shell  is  a function  of  the  L 
value.  ) The  reason  for  the  specification  of  the  50°  E.  drift  shell  will  be 

i 

made  clear  later.  We  can  further  divide  the  data  set  of  Fig.  5,  considering 

( 

only  those  points  which  are  defined  in  a)  and  b)  above.  We  then  find  that 
of  the  35  observation  points  prior  to  Day  262,  27  (or  77%)  exhibited  precipi- 
tation spikes.  Of  the  remaining  22  points  after  Day  262,  three  (or  14%) 
had  spikes.  Since  the  OVl-19  orbit  covered  only  a narrow  range  of  local 
time  (it  was  nearly  sun-synchronous),  possible  correlations  of  the  points 
of  observation  of  spike  events  with  local  time  could  not  be  ascertained. 

There  was  no  correlation  with  magnetic  activity. 

For  the  bigger  events,  sufficient  intensity  was  available  for  accurate 
determination  of  the  local  pitch-angle  distribution.  The  instrumentation 
contained  count-rate  meters  which  have  sufficiently  long  time-constants  at 
low  count-rates  that  accurate  determination  of  local  loss-cone  angles  was 
not  possible  for  most  of  the  events.  A peak  count- rate  of  at  least  200  counts/ 
sec  was  required.  For  these  larger  events,  one  can  construct  an  envelope 
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Figure  5 Scatter  plot  of  precipitation  events  (and  lack  thereof)  observed 

by  the  OVl-19  as  a function  of  East  Longitude  and  time.  Circles 
indicate  satellite  passes  in  which  the  altitude  was  below  the 
drift-loss-cone  defined  by  the  magnetic  field  geometry  at  50° 

E.  L.  An  abrupt  change  in  the  frequency  of  occurrence  of  events 
is  evident  about  Day  260. 


of  the  maxima  in  a spike  at  a given  energy.  Then,  using  L as  a parameter, 
one  can  normalize  all  the  data  points  in  the  event  at  that  energy  to  the  count- 
rate  at  the  central  L value,  e.  g. 


N 1 • j = N , , (exp[ ; L - L /k] } 

norm.alized  observed'-  center  ^ 


■•vhere  L is  the  instantaneous  L value  and  k is  derived  from  the  envelope. 

In  this  manner,  a fairly  accurate,  to  about  2 degrees,  determination  of 
the  local  loss-cone  angle  could  be  made.  The  procedure  also  involved 
obtaining  the  convolution  of  the  instrument  response  function  with  an 
assumed  loss-cone  angle.  Figure  6 gives  a typical  result  for  an  assumed 
loss-cone  angle  of  60°  convoluted  with  the  actual  angular  response  function 
for  the  channel  from  which  the  data  points  of  Fig.  6 were  obtained. 

The  data  points  are  normalized  data  points.  Note  that  the  data  must  be 
almost  two  orders  of  magnitude  above  background  before  the  time-constant 
effect  in  the  rate-meters  becomes  negligible.  For  the  data  of  Fig.  6,  the 
maximum  normalization  factor  was  about  6. 

Having  obtained  a reasonably  accurate  measurement  of  the  local 
loss-cone  angle,  the  equatorial  pitch-angle  of  the  observed  loss-cone  is 
easily  determined  using  the  local  B value  and  assuming  adiabatic  motion. 
This  angle  then  can  be  used  to  determine  the  location  of  the  origin  of  the 
spike  (see  Luhmann  and  Vampola,  1977  for  details).  If  one  assumes  that 
the  size  of  the  loss-cone  is  determined  by  the  size  of  the  local  bounce  loss- 
cone  at  the  point  of  origin,  one  can  trace  back  each  event  along  its  drift 
shell  to  the  east  longitude  at  which  the  instrumentally  determined  equatorial 
pitch-angle  of  the  loss  cone  coincides  with  the  local  bounce  loss  cone 
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Figure  6 A comparison  of  pitch-angle  data  from  one  event  with  the  instru 
mental  response  to  a hypothetical  distribution  which  is  isotropic 
between  60°  and  90°  and  zero  elsewhere.  The  data  points  have 
been  normalized.  The  "looping"  which  is  apparent  in  the  data 
at  the  low  fluxes  is  due  to  the  relatively  long  time-constants  of 
the  count-rate  meter  used  in  this  instrument. 


equatorial  pitch-angle.  We  cannot  show  that  the  pitch-angle  diffusion 
experienced  by  the  particle  after  its  initial  perturbation  is  negligible,  but 
we  believe  the  data  indicates  that  this  is  so. 

Figure  7 presents  the  data  of  Figure  3 as  an  "L"  vs  "Equatorial  Loss 
Cone  ^ngle"  plot.  For  this  display,  the  procedure  outlined  above  was  used 
for  each  of  the  energies.  For  the  139  keV.  312  keV.  and  380  keV  channels, 
a local  loss-cone  angle  of  90°  was  assumed,  since  the  pitch-angle  convolu- 
tion procedure  results  were  consistent  with  this  value,  although  for  the 
312  keV  channel  87°  or  88°  was  also  consistent.  The  other  two  channels 
had  local  loss-cones  of  84°  and  85°.  These  are  significantly  larger  local 
loss-cones  than  were  generally  observed  in  the  data  because  of  the  relatively 
low  altitude  of  the  observation.  For  most  of  the  events,  local  loss-cones 
of  50°  to  75°  were  observed.  However,  the  error  in  determination  of  the 
equatorial  loss-cone  angle,  is  considerably  less  for  a large  local  loss-cone 
angle.  For  instance,  at  89°.  a 1°  error  in  determination  of  the  local  loss- 
cone  angle  results  in  an  error  of  only  0.01°  in  the  equatorial  loss-cone  angle 
while  at  59°  local  loss-cone  angle,  a 1°  error  results  in  an  equatorial  loss- 
cone  angle  error  of  almost  0.  5°.  A 0.  5°  equatorial  angle  variation  translates 
to  a 5°  east  longitude  variation  for  values  around  60°  east  longitude,  as 
shown  by  Figure  4.  For  higher  altitude  data,  the  error  would  not  be  as 
extreme  as  this.  Figure  7 also  includes  a curve  showing  the  equatorial 
pitch-angle  of  particles  which  mirror  at  100  km  in  the  south  and  pass  through 
the  equatorial  region  at  60°  east  longitude.  Typical  mirror  longitudes  for 
these  particles  on  these  L-shells  are  66°  to  68°  in  the  south  and  59  m the 
north. 
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Figure  7 A plot  of  the  equatorial  pitch-angle  of  the  observed  local  loss- 
cone  vs.  L for  the  data  of  Figure  3.  The  solid  curve  gives  the 
equatorial  pitch-angle  as  a function  of  L for  particles  which  cross 
the  magnetic  equator  at  60°  East  Longitude  and  mirror  at  100 
km  in  the  southern  hemisphere.  Particles  mirroring  at  100  krr.  in  the 
northern  hemisphere  have  smaller  pitch-angles. 


Figure  8 shows  the  192  keV  data  for  sixteen  events  mapped  back 
from  their  points  of  observation  to  their  points  of  origin  as  defined  above. 
All  but  one  of  them  map  back  to  the  region  54°  to  62°  E.  with  half  of  those 
being  between  54°  and  56°.  An  error  in  determinatio  n of  the  local  loss 
cone  angle  of  1°  leads  to  an  error  of  about  0.  8°  in  east  longitude  for  these 
data. 

For  this  mapping,  we  assumed  the  loss-cone  angle  was  produced  by 
a minimum  altitude  point  in  the  southern  hemisphere.  vVe  then  determined 
the  east  longitude  of  the  field-line  which  had  the  appropriate  bounce-loss 
cone  equatorial  pitch-angle,  and  mapped  back  to  the  equatorial  crossing 
along  the  field  line.  The  field  model  used  was  the  POGO  (8/691  m.odel  (Cain 
and  Sweeney,  19T0). 

The  data  of  Figure  7 indicate  that  particles  of  various  energies  are 
being  precipitated  at  the  sam.e  longitude  but  not  at  the  sam.e  L shell.  The 
data  of  Figures  2 and  3 indicate  both  regularity  and  monotonicity  in  the 
distribution  of  the  various  energy  peaks  as  a function  of  L.  To  investigate 
this  possible  correlation,  the  events  used  in  Figure  8 to  map  back  in 
longitude  were  also  plotted  on  an  ENERGY  vs.  L diagram,  shown  as 
Figure  9.  For  this  plot,  the  L value  of  the  peak  intensity  for  each  energy 
channel  was  used,  with  an  accuracy  of  about  = .01  L.  The  energy  assigned 
to  each  channel  was  the  centroid  of  the  energy  response  for  a flat  spectrum 
( which  should  be  fairly  well  approximated  at  the  peak  of  each  precipitation 
spike,  as  shown  by  Figure  1).  Calculations  of  the  expected  resonant  energy 
at  the  equator  as  a function  of  L for  a 16.  2 kHz  wave  (non-ducted)  for 
two  typical  cold  electron  densities  are  also  shown.  These  calculations 
were  done  by  the  m.ethod  of  ray- tracing  using  a diffusive  equilibrium 
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Figure  8 A mapping  of  observed  precipitation  events  back  to  their  pre- 
sumed origin  using  the  method  outlined  in  the  text.  These  are 
all  OVl-19  data  from  the  192  keV  channel.  The  large  variation  in 
equatorial  loss-cone  angle  mapping  back  to  approximately  the 
same  East  Longitude  is  due  to  a relatively  large  variation  in  L 
from  event  to  event  for  the  peak  in  precipitation. 
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Figure  9 L dependence  of  the  location  of  the  peak  at  each  energy  for  \ 

the  events  of  Figure  8.  The  dashed  lines  are  the  expected  J 

resonant  energies  at  the  equator  for  a 16.  2 kHz  wave  for  two  i 

typical  diffusive -equilibrium  electron  density  models.  Nq  j 

refers  to  the  base  level  electron  density  at  1000  km.  I 
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model  for  the  magnetospheric  electron  density  model  of  Angerami  and 
Thomas  (1964).  The  model  parameters  consisted  of  a constant  base 
electron  density  at  1000  km  and  ion  mixture  of  90%/10%  at  a tempera- 

ture of  1600  K.  At  1000  km  for  an  L of  1.97  during  the  spring  of  1965, 

Brace  et  al.  (1967)  reported  variations  from  1.3  x 10*  cm”  at  a local 

j.  - 3 

time  of  1300  hr s to  2.5  x 1 0 * cm  at  0 100  hr  s LT.  These  numbers  are 
quite  similar  to  the  values  used  in  Figure  9. 


-Z4  - 


DISCUSSION 


The  data  of  Figures  1-3  indicate  a narrow-band  resonance  is  in- 
volved in  the  precipitation  of  these  electrons  in  the  inner  zone.  While 
the  possibility  of  such  a narrow  resonance  being  caused  by  a natural 
instability  cannot  be  ruled  out  on  fundamental  grounds,  no  such  narrow 
resonance  has  been  observed  in  wave  data  from  satellites.  (We  are 
not  considering  the  various  cold  proton  and  electron  gyroresonances  as 
possibilities  here;  nor  are  we  going  to  consider  anomalous  amplification 
of  VLF  signals  as  being  a natural  source.  ) 

If  one  makes  the  assumption  that  the  precipitations  are  due  to 
artificially  generated  VLF  waves,  the  first  task  is  to  identify  the  source 
of  those  waves.  The  mapping  procedure  which  resulted  in  Figure  8 
indicates  that  the  source  of  all  the  precipitations  examined  was  west  of 
80°  E.  L.  and,  for  all  but  one  event,  was  west  of  61°  E.  L.  Figure  4 
indicates  that  by  the  time  the  electron  population  has  drifted  from  330°  E. 
to  60°  E.  , the  L = 1.75  bounce  loss-cone  angle  decreases  from  about 
28°  to  22°  (equatorial  pitch-angle).  The  highly  asymmetric  pitch- angle 
distribution  created  by  the  South  Atlantic  Anomaly  seems  to  be  a prime 
candidate  for  a Kennel-Petschek  (1966)  instability.  The  empty  portion  of 
the  drift  loss  cone  would  serve  as  an  excellent  background  against  which 
to  see  such  an  instability  precipitate  a portion  of  the  asymmetric  distri- 
bution. However,  we  can  rule  that  out  as  a source  for  these  precipitations 
on  a number  of  grounds. 

First  of  all,  such  a process  would  not  be  energy  selective,  since 
the  asymmetric  distribution  occurs  at  all  energies.  Figure  1 shows  our 
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I process  ^ energy  selective.  Second,  the  process  would  not  be  highly 

selective  in  L.  Figures  2 and  3 again  show  that  ours  is.  TViird,  the 
data  of  Figure  5 show  an  abrupt  change  in  the  frequency  of  occurrence 
of  these  events  after  Day  260.  This  is  also  incompatible  v^ith  a K-P 
instability  as  the  source.  Finally,  there  was  no  correlation  bepveen 

f these  events  and  magnetic  activity.  Since  the  low  energy  particle  pcpu- 

i I 

lation  at  these  L values  is  enhanced  during  ancl  following  major  magnetic 

' storms  ( Bostrom  et  al.,  1970;  Vampola,  1972),  the  asvm.m.etry  east  of 

the  anomaly  is  greater  follo'ving  such  storms  than  before.  Large  storms 
occurred  around  Day  135  and  Day  270,  but  these  had  no  effect  on  the 
frequency  of  observance  of  events  or  on  the  am.ount  of  precipitation  occur- 
ring in  a single  event. 

The  source  of  these  precipitation  events  rotates  with  the  earth,  as 
shown  by  the  fact  that  15  of  the  16  events  mapped  in  Figure  8 originated 
between  54°  and  62°  E.  Figure  7,  for  instance,  shows  that  all  energies, 
although  at  different  L values  and  with  different  equ?torial  loss  cone 
angles,  last  interacted  strongly  w'ith  the  atmosphere  at  about  60°  E.  Since 
the  equatorial  loss-cone  angle  is  decreasing  with  eastward  drift  at  that 
location,  we  can  assume  that  all  energies  were  precipitated  there.  Of 
course,  particles  could  also  have  been  precipitated  at  longitudes  farther 
west,  but  those  particles  would  have  a larger  loss-cone  angle.  The  fact 
I that  events  observed  half-way  around  the  world  to  the  east  from  where  they 

last  interacted  with  the  atmosphere  show  similar  loss-cone  angles  indi- 
cates that  the  pitch-angle  diffusion  due  to  non-atmospheric  causes  is  very 
weak.  Otherv5,'ise,  there  would  be  a tendency  for  events  observed  farther 
east  to  have  smaller  loss-cones.  The  slightly  smaller  loss-cones  observed 
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at  lower  energy  in  Figure  7 may  be  due  either  to  slow  pitch-angle 
diffusion  which  is  more  effective  at  low  energy  or  may  be  indicative 
that  the  precipitation  mechanism  continues  to  operate  farther  to  the 
east  for  lower  energy  particles.  A precipitation  cause  which  rotates 
with  the  earth  is  strongly  suggestive  of  a ground-based  VLF  transmitter. 
The  radiation  pattern  from  the  transmitter  diminishes  in  intensity  with 
distance  and  hence  might  be  effective  only  on  the  lowest  energy  electrons 
at  greater  distances. 

In  the  neighborhood  of  60°  E.  , there  were  at  least  two  VLF  sta- 
tions: a relatively  low  power  transmitter  at  Algazy,  USSR,  at  77°  11'  E. 
43°  50'  N.  that  operated  at  16.  1 kHz  (Aksenov  et  al. , 1970),  and  a high 
power  (315  kw)  transmitter  at  Gorky,  USSR,  at  44°  E.  , 46.2°  N.  , that 
has  operated  as  UMS  on  17.  1 kHz  and  as  RPS  at  16.2  kHz  (Watt,  1967; 
Amon  and  Dowden,  1976). 

The  plot  of  the  energy  dependence  of  the  location  of  the  peaks  as 
a function  of  L shown  in  Figure  9 indicates  that  a single  frequency  with  a 
single  electron  density  for  the  unducted  calculations  gives  reasonable 
agreement.  The  slope  is  similar  to  that  reported  by  Imhof  et  al.  (1974a) 
for  L-dependent  peaks  in  1971-1972.  Frequencies  of  8 to  1 5 kHz  were 
indicated  as  being  compatible  with  those  observations.  For  those  data, 
Imhof  et  al.  concluded  that  the  peaked  spectra  precipitations  occurred 
preferentially  in  the  early  dawn  hours.  In  our  data,  we  consider  the 
variations  as  being  due  to  a local-time  variation  in  the  electron  densitv' 
with  a constant  perturbing  frequency  rather  than  a variable  perturbing 
frequency.  Detailed  wave-tracing  studies  on  this  data  set  are  in  progres 


The  sudden  change  in  the  frequency  of  occurrence  of  these  precipitation 
spikes  which  occurred  after  Day  260  is  currently  unexplained.  If  a trans- 
mitter at  Gorky  is  assumed,  precipitation  spikes  should  be  observable 
whenever  the  satellite  is  in  the  drift  loss-cone  above  an  altitude  determined 
by  the  bounce  loss-cone  at  the  point  farthest  to  the  east  that  the  precipita- 
tion mechanism  is  effective.  Since  most  of  the  events  m.apped  back  to 
55°-60°,  a conservative  limit  of  50°  was  set  for  the  determination  of  the 
presence  or  absence  of  an  event.  Previous  studies  of  triggered  emissions 
indicate  that  a CW  signal  of  about  125  msec  is  required  to  trigger  the 
instability  that  leads  to  wave  growth  (and  presumably  particle  pitch-angle 
scattering).  If  a single  dash  in  a transmission  from  Gorky  can  produce 
precipitation  to  as  far  as  60°  E.  L. , a region  of  tens  of  degrees  in  longitude 
presumably  is  filled  with  drift  loss-cone  electrons.  Since  the  drift  rate 
of  these  electrons  is  of  the  order  of  2°  per  minute,  it  would  take  on  the 
order  of  ten  minutes  for  the  effects  of  a single  dash  to  drift  through  the 
plane  of  the  satellite's  orbit.  Hence,  relatively  infrequent  communica- 
tions transmissions  would  be  sufficient  to  produce  effectively  steady- 
state  precipitation  effects.  However,  it  is  known  that  the  transmissions 
from  Gorky  were  not  always  GW.  Frequency  Shift  Keying,  which  was  also 
used,  probably  would  not  produce  precipitation.  The  duration  of  transmission 
of  a given  frequency  in  FSK  is  probably  too  short  to  produce  triggered 
emissions,  ^n  additional  piece  of  evidence  of  the  efficacy  of  producing 
particle  precipitation  by  FSK  is  given  indirectly  by  BuUough  et  al.  (1976^ 
in  which  GBR  signals  at  16  kHz  received  aboard  Ariel  III  were  character- 
istically an  order  of  magnitude  stronger  in  GW  than  in  FSK.  Hence,  the 
change  in  frequency  of  observation  of  precipitation  events  which  occurred 
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during  1969  may  reflect  a change  in  the  operating  mode  of  the  transmitter. 
A complete  study  of  this  effect  would  require  a detailed  operating  log  from 
the  transmitter.  Such  a log  would  permit  quantitative  analysis  of  the 
effe  :ts,  including  the  drift  periods  and  the  local  pitch-angle.  It  would 
also  permit  identification  of  the  area  over  which  the  transmissions  were 
able  to  produce  precipitation. 

In  addition  to  the  events  included  in  this  study,  other  precipitation 
spikes  with  different  characteristics  were  also  observed.  Several  were 
"multiple  peak"  events,  with  two  or  more  peaks,  at  different  L,  for  each 
energy.  Possible  explanations  are:  a)  multiple  transmitting  frequencies, 
perhaps  as  much  as  10  minutes  apart  in  time;  b)  transmitters  at  different 
sites  (and  therefore  at  different  local  times,  providing  different  base 
electron  densities);  c)  a non-uniform  ionosphere.  Also,  a single  "ducted" 
event  was  observed  in  which  all  channels  up  to  and  including  the  531  keV 
channel  were  observing  precipitation  at  the  same  L shell.  The  observed 
loss-cone  coincided  with  the  local  loss-cone.  Hence,  either  the  duct  was 
local  (at  288°  E.  L.  ) or  could  have  been  located  anywhere  east  of  about 
40°  E.  L.  All  these  events  are  being  analyzed. 

The  practical  and  geophysical  implications  of  artifically  induced 
precipitation  of  radiation  belt  electrons  are  far-reaching.  Bullough  et  al. 
(1976)  suggested  that  harmonics  of  the  extensive  50-60  Hz  power  grids 
in  North  American  and  Europe  are  effective  in  precipitating  electrons 
and  may  be  responsible,  in  part,  for  the  existence  of  the  slot  at  2 < L < 3. 
Imhof  et  al.  (1974b)  observed  what  they  called  "resonance  precipitation" 
of  electrons  in  the  slot  region  which  was  consistent  with  plasmaspheric 
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hiss  observed  in  the  slot.  The  hiss  may  be  the  result  of  amplification  of 
high  order  harmonics  of  power  grids.  If  other  man-made  sources  of 
ELF/VLF  are  effective  in  the  precipitation  of  electrons,  and  there  are  a 
number  of  very  powerful  transmitters  around  the  world  used  for  naviga- 
tion and  communications,  the  possibility  exists  that  the  gross  morphology 
of  the  radiation  belts  may  be  the  result  of  man's  transmissions.  Perhaps 
weak  pitch- angle  diffusion  observed  in  the  outer  belt  (e.  g.  , Williams  et  al.  , 
1968)  is  entirely  due  to  such  emissions.  A significant  and  unambiguous 
confirmation  that  this  mechanism  is  operative  could  be  made,  of  course, 
by  turning  off  all  sources  of  man-made  ELF/VLF/LF,  including  power 
grids.  A realistic  method  would  be  to  produce  time- averaged  intensity 
plots  of  particles  in  the  bounce  loss-cone  as  a function  of  longitude  and 
latitude.  If  only  a small  number  of  discrete  sources  of  precipitation  exist, 
they  would  be  readily  observable.  The  effect  of  the  power  grids  over 
North  America  should  be  discernable  if  sufficient  resolution  in  both  energy 
and  latitude  is  available. 
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SUMM-ARY 


Highly  monoenergetic  spikes  of  energetic  electrons  were  observed 
in  the  drift  loss-cone  at  L-values  of  1.  6 to  1. 85  during  the  1968-  1970 
period  by  magnetic  spectrometers  on  the  OVl-14  and  OVl-19  satellites. 

Using  the  locally  observed  loss-cone  angles,  these  spikes  were  determined 
to  have  last  interacted  with  the  atmosphere  in  the  region  of  55°  to  62°  East 
Longitude,  indicating  that  the  cause  of  the  precipitation  was  corotating 
with  the  earth.  The  energy  dependence  of  the  position  in  L of  the  precipi- 
tation was  shown  to  agree  with  the  predicted  interaction  region  for  unducted 
waves  with  a frequency  similar  to  the  Russian  communications  VLF  trans- 
mitter which  was  located  near  the  observed  source  of  the  precipitations. 

-Although  the  geomagnetic  topography  over  eastern  Europe  makes  this  an 
ideal  location  for  observing  effects  of  VLF  transmissions  on  the  trapped 
particle  population,  since  in  this  region  the  South  -Atlantic  -Anomaly  has 
produced  a large  empty  loss-cone,  it  should  be  possible  to  see  similar 
effects  from  other  powerful  VLF  transmitters  by  carefully  studying  the 
local  bounce  loss-cone  on  a worldwide  basis.  The  identification  of  direct  1 

precipitation  of  electrons  by  ground-based  VLF  stations  adds  further 
support  to  the  idea  advanced  by  Bullough  et  al.  (1976)  that  significant  fea- 
tures in  the  electron  belts  may  be  due  to  man-made  transmissions. 
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THE  IVAN  A.  GETTING  LABORATORIES 

The  Laboratory  Operations  of  The  Aerospace  Corporation  is  conducting 
experimental  and  theoretical  investigations  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems.  Ver- 
satility and  flexibility  have  been  developed  to  a high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  scientific  devel- 
opntents  is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are; 

Aerophysics  Laboratory;  Launch  and  reentry  aerodynamics,  heat  trans- 
fer, reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high-power  gas  lasers. 

Chemistry  and  Physics  Laboratory:  Atmospheric  reactions  and  atmos- 
pheric optics,  chemical  reactions  in  polluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo- 
sensitive materials  and  sensors,  high  precision  laser  ranging,  and  the  appli- 
cation of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory;  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plasma  electromagnetics;  quantum  electronics, 
lasers,  and  electro-optics;  communication  sciences,  applied  electronics,  semi- 
conducting, superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging;  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 

Materials  Sciences  Laboratory;  Development  of  new  materials;  metal 
matrix  composites  and  new  forms  of  carbon:  test  and  evaluation  of  graphite 
and  ceramics  in  reentry;  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment;  application  of  fracture  mechanics  to  stress  cor- 
rosion and  fatigue-induced  fractures  in  structural  metals. 

Space  Sciences  Laboratory;  Atmospheric  and  ionospheric  physics,  radia- 
tion from  the  atmosphere,  density  and  composition  of  the  atmosphere,  aurorae 
and  airglow;  magnetospheric  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere;  solar  physics,  studies  of  solar  magnetic 
fields;  space  astronomy,  x-ray  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere;  the  effects  of  optical,  electromagnetic,  and  particulate  radia- 
tions in  space  on  space  systems. 
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